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Abstract 

Wc report the results of a study of neutral B decays to the D^tt^tt~ and D*^tt^tt~ final states 
using complete D^*^'^ reconstruction. The contributions from two-body B D**7r with narrow 
{jq = 3/2) D** states and B D^*^ p, D^*^ f2, D^*^a decays have been determined. All results are 
preliminary, and are based on a large data sample collected in the Belle experiment at the KEKB 
e+e~ collider. 

PACS numbers: 13.25.Hw, 14.40Lb, 14.40.Nd 
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The decays of B meson to Dir and D*7r final states are two of its dominant hadronic 
decay modes and have been measured quite well 0. In this paper we study the production 
of excited states of D-mesons, collectively referred to as D**'s, that are P-wave excitations 
of quark-antiquark systems containing one charmed and one light {u, d) quark. The results 
provide tests of Heavy Quark Effective Theory (HQET) and QCD sum rules. FigureHshows 
the spectrum of D-meson excitations. In the heavy quark limit, the heavy quark spin Sc 
decouples from the other degrees of freedom and the total angular momentum of the light 
quark jq = L+Sg is a good quantum number. There are four P-wave states with the following 
spin-parity and light quark angular momenta: 0+(jg = 1/2), l^{iq = 1/2), l^{jq = 3/2) 
and 2~^{jq = 3/2), which are usually labeled as D^, D[, Di and D2, respectively. 




FIG. 1: Spectrum of D-meson excitations. Lines show possible one pion transitions. 

The two jn = 3/2 states are narrow with widths of about 20-40 MeV and were previously 
observed 0, 0, B la, IE B i j l , | R , 10- '111 1^ . The measured values of their masses agree with 



model predictions [l3|, |l5|, |l6| . The remaining jq = 1/2 states decay via S- waves and are 
expected to be quite broad. 

The B D^*^7m decay provides a possibility to study D** production. The fixed spin of 
the initial state makes it possible to perform an angular analysis of the decay products and 
to separate final states with different quantum numbers. 

We have acquired a large sample of B decays using the Belle detector, which has good 
resolution and particle ID. These data can be used to contribute to an understanding of 
these problems. 

Earlier neutral D** production in charged -B-decays has been studied at Belle AH 
four D** states have been observed and the production rate of the broad (j = l/2)-states 
appears to be comparable to that of the narrow (j = 3/2)-states. This work describes a 
similar analysis of the decay B^ D**^tt~ . 
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In the case of neutral B decay the D^*^7m final state contains two pions of opposite sign 
that can form several bound states (such as p, /o, /2), which should also be taken into 
account. The presence of vrvr bound states complicates the analysis but can give valuable 
information about the constants and mechanism of these decays. 



THE BELLE DETECTOR 



The Belle detector [18[ is a large-solid-angle magnetic spectrometer that consists of a 
three-layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC) for charged 
particle tracking and specific ionization measurement (dE/dx), an array of aerogel threshold 
Cerenkov counters (ACC), time-of-fiight scintillation counters (TOF), and an array of 8736 
CsI(Tl) crystals for electromagnetic calorimetry (ECL) located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron fiux return located outside the 
coil is instrumented to detect mesons and identify muons (KLM). We use a GEANT- 
based Monte Carlo (MC) simulation to model the response of the detector and to determine 
its acceptance |l9l |. 

Separation of kaons and pions is accomplished by combining the responses of the ACC 
and the TOF with dE/dx measurements in the CDC to form a likelihood C{h) where h = (vr) 
or (K). Charged particles are identified as pions or kaons using the likelihood ratio (TZ)'- 

At large momenta (>2.5 GeV/c) only the ACC and dE/dx are used since the TOF provides 
no significant separation of kaons and pions. Electron identification is based on a combina- 
tion of dE/dx measurements, ACC photoelectron yields and the position, shape and total 
energy deposition (E/p) of the shower detected in the ECL. A more detailed description of 
the Belle particle identification can be found in Ref. [2o| . 



EVENT SELECTION 

A data sample of 140 fb^^ (152 million BB events) collected at the T{4S) resonance with 
the Belle detector is used. Candidate B^ — > D^7r~^7r~ and B^ — >■ D*^7r~^'ir~ events as well 
as charge conjugate combinations are selected. The and D*'^ mesons are reconstructed 
in —>■ K~7[~^ and D*^ D^n^ modes, respectively. from D*^ decay is detected in 
channels K~tt^ and K'h^tt^h' . The signal-to-background ratios for other D 

decay modes are found to be much lower and they are not used in this analysis. 

Charged tracks are selected with requirements based on the average hit residuals and 
impact parameters relative to the interaction point. We also require that the polar angle of 
each track be within the angular range of 17° — 150° and that the transverse track momentum 
be greater than 50 MeV/c for kaons and 25 MeV/c for pions. 

Charged kaon candidates are selected with the requirement TZ{K) > 0.6. This has an 
efficiency of 90% for kaons and a pion misidentification probability of 10%. For pions the 
requirement 7?.(7r) > 0.2 is used. All tracks that are positively identified as electrons are 
rejected. 

mesons are reconstructed from K'tt^ combinations with invariant mass within 
12 MeV/c^ of the nominal mass, which corresponds to about 2.5 aKn- We reconstruct 1?*° 
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mesons from the Dvr^ combinations with a mass difference of Mp^o — M£)0 within 2.5 MeV/c^ 
of its nominal value. 

Candidate events are identified by their center of mass (cm.) energy difference, AE = 



iJ2iEi) — Ey,, and beam-constrained mass, Mbc = y E^ — {J^iViYi where E\, = \/s/2 is the 
beam energy in the T(4S') cm. frame, and pi and Ei are the cm. three- momenta and energies 
of the B meson candidate decay products. We select events with Mbc > 5.25 GeV/c^ and 
\AE\ < 0.10 GeV. 

To suppress the large continuum background (e+e^ qq, where q = u,d, s, c), topological 
variables are used. Since the produced B mesons are almost at rest in the cm. frame, the 
angles of the decay products of the two B mesons are uncorrelated and the tracks tend to 
be isotropic while continuum qq events tend to have a two-jet structure. We use the angle 
between the thrust axis of the B candidate and that of the rest of the event (Othrust) to 
discriminate between these two cases. The distribution of | cos 6thrust I is strongly peaked 
near | cosGthrustI = 1 for qq events and is nearly flat for T(4S') BB events. We require 
I cosGthrustI < 0.8, which eliminates about 83% of the continuum background while retaining 
about 80% of signal events. 

There are events for which two or more combinations pass all the selection criteria. 
According to a MC simulation, this occurs primarily because of the misreconstruction of a 
low momentum pion from the D** D^*^n decay. To avoid multiple entries, the combination 
that has the minimum difference of Z coordinates at the interaction point, —Z-j^^ly of the 



tracks corresponding to the pions from B — > D**tti and D** — >■ D^*'^Tr2 decays is selected 21 



This selection suppresses the combinations that include pions from Ks decays. In the case 
of multiple D combinations, the one with invariant mass closest to the nominal value is 
selected. 



B° D^TT+TT~ ANALYSIS. 

The final state of the B^ —>■ D^-k^-k^ decay together with three-body and quasi-two-body 
events includes the two-body decay B'^ — *• D*^7r~ followed by the decay D*^ D^tt^ . 
Using the mass difference of M/)^ — M/j we subdivide the total sample in two: events with 
IMdtt - Md- 0.1455| < 0.03 (GeV/c^)(~ Ga) (denoted further as sample (2)) correspond 
to D*7r production and the rest of the events Dt^tt form sample (1). 

The Mbc and AE distributions for B^ — > D^tt^tt' events are shown in Fig. |21 The 
distributions are plotted for events that satisfy the selection criteria for the other variable: 
i.e., \AE\ < 25 MeV and |Mbc - Mb\ < 5 MeV/c^ for the Mbc and the AE histograms, 
respectively. A clear signal is evident in both distributions. The signal yield is obtained 
by fitting the AE distribution to the sum of two Gaussians with the same mean for the 
signal and a linear function for background. The widths and the relative normalization of 
the two Gaussians are fixed at values obtained from the MC simulation while the signal 
normalization as well as the constant term and slope of the background linear function are 
treated as free parameters. 

The signal yields are 1128 ± 51 and 1521 ± 40 for samples (1) and (2), respectively. A 
detection efficiency of (16.8 ± 0.4)% and (18.2 ± 0.4)% is determined from a MC simulation 
that uses a Dalitz plot distribution that is generated according to the model described in the 
next section. Taking into account the branching fraction i3(i^° K^n^) = (3.80 ± 0.09)% 
and B{D*^ D^tt^) = (67.7 ± 0.5)% |]J, we obtain the following value for the branching 
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data AE(GeV) data AE (GeV) 

FIG. 2: Mbc and AE distributions for S° vr+TT events. (a),(b) show Mbc distributions. 

(c),(d) show AE distributions. (a),(c) ((b), (d)) distributions are plotted for sample (1)((2)), re- 
spectively. 

fractions: 

B{B° D^TT+TT-) = (1.07 ± 0.06 ± 0.10) x 10'^, 

and 

B{B° D*+7T-) = (2.30 ± 0.06 ± 0.19) x 10'^, 

where the first error is statistical and second error is systematic. The main contributions 
to the systematic error are hsted in Table IH The uncertainty of the background shape was 
estimated by adding higher order polynomial terms to the approximating function. 
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MC 
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2.4% 


2.5% 


Total 


9.2% 


8.1% 



TABLE I: The systematic uncertainties for — > D^tt^-k . 

The values of B{B'^ D^tt^tt') and B{B^ D*+7r^) are in agreement with a previous 
result from Belle: ^ D^tt+tt-) = (8.0 ± 1.6) x 10"^ 22] and a CLEO result: ^ 

D*+n-) = (2.81 ± 0.25) x 10"^ ^ and are more precise. ' 



i? — > Dtttt Dalitz plot analysis 

For a three-body decay of a spin zero particle, two variables are required to describe 
the decay kinematics; we use the D^vr"*" and vr+vr" invariant masses squared, and qf, 
respectively. 

To analyze the dynamics of S — > Dtttt decays, events with AE and Mbc within the 
signal region {{AE + ^(Mbc - MB))/cr^Ey + ((Mbc - MB)/aMbJ^ < 4 are selected. The 
parameters ctae = llMeV/c^, 0"^//^^ = 2.7 MeV/c^, k = 0.9 have been obtained from a fit to 
experimental data; the coefficient k takes into account a correlation between Mbc and AE. 

To model the contribution and shape of the background, we use events from the AE side- 
bands, which are defined as: {{AE ±65 MeV + k{M^^ - Mb)) / a aeY + {{M^^- Mb)/ a m,,,)'^ < 
4. Figure El shows the signal and sidebands regions in the (Mbc — ^E)) plane. 

The Dtt and tttt mass distributions for the signal and sideband events are shown in 
Fig. m In the Dti mass distribution we can clearly see the narrow peak of the The tttt 
distribution has a signal of the p meson as well as a structure at 1.2 — 1.3 GeV/c^ that can 
be due to /o(1370) or /2(1270) contributions. 

The distributions of events in the M|,^ versus M^^ Dalitz plot for the signal and sideband 
regions are shown in Fig. El The Dalitz plot boundary is determined by the decay kinematics 
and the masses of the daughter particles. In order to have the same Dalitz plot boundary for 
events in both signal and sideband regions, mass-constrained fits of Ktt to M/) and Dtttt to 
Mb are performed. The mass-constrained fits also reduce smearing from detector resolution. 

To extract the amplitudes and phases of different intermediate states, an unbinned fit 
to the Dalitz plot is performed using the method described in Ref. The event density 
function in the Dalitz plot is the sum of the signal and background. 

The background distribution and normalization are obtained from the AE sideband anal- 
ysis. Since the Dtt mass distributions for the upper and lower halves of the AE sideband 
have similar shapes, we can expect similar background behavior for the signal and sideband 
regions. The background Dalitz plot has neither a resonant structure nor non-trivial he- 
licity behavior and is combinatorial in its origin. The background shape is obtained from 
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" *5.25 5.26 5.27 5.28 5.29 

M,^ (GeV/c') 

FIG. 3: The experimental event distribution in the (Mbc — ^E)) plot. The ellipses show a position 
of the signal (1) and sideband regions (2), (3). 



an unbinned fit of the sideband distribution to a smooth two-dimensional function. The 
number of background events in the signal region is scaled according to the relative areas of 
the signal and the sideband regions. 

There is no generally accepted way to exactly describe a three-body amplitude. In this 
paper we represent the Dt^tt amplitude as the sum of Breit-Wigner contributions for different 
intermediate two-body states. Such an approach cannot be exact since it is neither analytic 
nor unitary and does not take describe completely possible final state interactions. Never- 
theless, the sum of Breit-Wigners describes the main features of the amplitude behavior and 
allows one to find and distinguish the contributions of various two-body intermediate states, 
their interference and the effective parameters of these states. We used the same approach 
in the analysis of charged B decays jl7l |. 

In the D^vr+vr^ final state a combination of the D'^-meson and a pion can form a vector 
meson a tensor meson or a scalar state -Dg^; the axial vector mesons Di and D'l' 
cannot decay to two pseudoscalars because of angular momentum and parity conservation. 
The region of the D^ti^ invariant mass corresponding to the D*^ is excluded from the 
fitting by requiring \M£,-k — M/)*! > 0.01 GeV/c^ but in B decay a virtual D*^ (referred 
to as Dl) can be produced off-shell with ^/q^ larger than the Z)°7r+ total mass and such a 
process will contribute to the amplitude. Another virtual hadron that can be produced in 
this combination is B*~ (referred to as B*): B —>■ B*tt and B* — »• Dir. For the mass of 
B*~ as well as the mass and width of D*~^, we used the PDG values 0; the widths of B*~ 
are calculated from the width of the D*~ in the HQET approach. In the tttt distribution 
we can see a p meson signal and some evidence of resonances around 1300 MeV, which we 
describe with the hypothesis of /2(1270) and /o(1370) mesons. As it is shown in Table ITTl 
the hypothesis of /2 has the best likelihood value. We also include the scalar resonance 
/o(600) with free mass and width. 

The contributions from the intermediate states listed above are included in the signal- 
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FIG. 4: Dtt - (a) and vrvr -(b) mass distribution. The points with error bars correspond to the 
signal box events, the hatched histogram — to the background obtained from sidebands. The open 
histogram shows the fit function after efficiency correction. 
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TABLE II: Comparison of the models with different resonances. 

event density (5'(g^, qf)) parameterization as a coherent sum of the corresponding amphtudes 
together with a possible constant amphtude (ops): 

S{q', ql) = \a,A\q\ ql) + a^e^^^ A\q\ ql) + a.e^^^ A\q\ ql) 
+ape''^'AP + aj^e'^^h+'t'p'^Afo + aj^e'^'^'h+'t'p) A^^ 
+aB*e"^^A\q\ql) + apse"^-^\\ (1) 

Each resonance is described by a relativistic Breit-Wigner with a q^ dependent width and 
angular dependence that corresponds to the spins of the intermediate and final state particles 
following the approach described in jT^J. We take into account hadronic transition form 
factors. The Blatt-Weisskopf parameterization jl^ with a hadron scale r=1.6 (GeV/c)~^ is 
used. For the virtual mesons and 5* that are produced beyond the peak region, another 
form factor parameterization is used: 

Fab{y>) = e-^^P-P"); (2) 
this provides stronger suppression of the Breit-Wigner far from the resonance region. 
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FIG. 5: The Dalitz plot of a)signal events; b)sideband events. 



The detector resolution for the invariant mass of the D-Kiini) combination is about 2.5 
(3.5) MeV which is much smaller than the narrowest peak width 30-40 MeV. 

The masses and widths of the tttt resonances (except /o(600)) are fixed at the PDG 



values. The mass and width of the broad resonance in (Dn): M 



2308MeV/c2, 1° 



276MeV/c^ have been taken from our measurement for the D**^ jl7| . 

The masses and widths of the D2 and the /o as well as the relative amplitudes and phases 
are free parameters of the fit. The parameters of the the /o are quite uncertain and it can 
also be regarded as a nonresonant S-wave structure. 

Table IIIII gives the results of the fit for different models. The contributions of different 
states are characterized by the branching fractions, which are defined as: 



Bn 



(3) 



where Ai{Q) is the corresponding amplitude, aj and (pi are the amplitude coefficients and 
phases obtained from the fit. The integration is performed over all the available phase 
space characterized by the multidimensional vector Q (for decay to 3 spinless particles 
dQ = dq'^dqf), and i is one of the intermediate states: D2, Dq, p, /2, /o, D*, B* or the 
constant term aps- 

In addition to the main minimum where the likelihood value is maximal Cq, there are 
several local minima with smaller likelihood C The value of (—2 In C) differs from (—2 In Cq) 
by 1.9-25. To search for local minima we perform 100 minimizations starting from different 
points, randomly distributed in the space of the minimized parameters. Local minima appear 
mainly as a result of different phases between the Dn and tttt structures. We consider the 
spread of the branching values as a model error. The central values are the parameters 
obtained for the main minimum. The spread of the relative phases is rather large. In some 
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0.44± 0.25 














l.OOib 0.21 












0.06± 0.05 




Mf„{GeV/c') 


0.658±0.062 


0.681±0.048 


0.649±0.056 


0.633±0.052 


0.676±0.057 


Tf„{GeV/c^) 


0.94± 0.22 


0.72± 0.14 


l.OOib 0.20 


0.88± 0.19 


0.78± 0.17 



TABLE III: The fit results for different sets of amplitudes. 



cases, the phases '"flip" by close to tt radians, so that the extraction of relative phases 
unreliable. 

The values of the resonance mass and width obtained from the fit are: 

M^.+ = (2459.5 ± 2.3 ± 0.7+^:^)MeV/c^ r^*+ = (48.9 ± 5.4 ± 4.2 ± 1.9)MeV. 

These parameters are consistent with the measurements from the FOCUS experiment: 
Md*o = (2467.6 ± 1.5 ± 0.8) MeV/c^, Yd* = 34.1 ± 6.5 ± 4.2) MeV 0. 

The product of the branching ratios for D2 production obtained in the analysis when the 
amplitudes are included from above is: 

B{B^ D*+7i-) X B{D;+ D%+) = (3.08 ± 0.33 ± Omtoll) x 10^^ 

where the indicated errors are the statistical, systematic and model error. 
The broad resonance branching fraction assuming 0+ quantum numbers is: 

B{B° D*+7i-) X B{D*+ D%+) = (0.60 ± 0.17 ± 0.16+°;^?) x 10"^ 

The helicity of the Dn distribution for different regions of is shown in Fig. [Altogether 
with the efficiency corrected fitting function. The histogram in the region of the D2 meson 
clearly indicates a D-wave. The distributions in the other region show reasonable agreement 
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of the fitting function and the data but the present statistics does not allow confirmation 
of the quantum numbers of the resonance. Table |n] shows that the likelihood changed 
significantly if we exclude this resonance but differs only slightly if the broad resonance is 
replaced with a vector or a tensor. We thus set an upper limit for the branching fraction of 
the scalar Dq-k^ production: 

^ D*+Ti-) X B{D*+ D\+) < 1.2 X 10"^ at 90% C.L.. 

The uncertainty of the background is one of the main sources of the systematic errors. 
It is estimated by comparing the fit results for the case when the background shape is 
taken separately from the lower or upper sideband in the AE distribution. The fit is also 
performed with more restrictive and loose cuts on AE, Mbc and AM^ that changes the 
signal-to-background ratio by more than a factor of 2. The obtained results are consistent 
with each other. The maximum difference is taken as an additional estimate of the systematic 
uncertainty. For the branching fractions, the systematic errors also include uncertainties on 
track reconstruction and PID efficiency, as well as the error in the K~7r~^ absolute 

branching fraction. 

The model uncertainties are estimated by comparing fit results for the case of different 
models and for the values of the parameter r of the transition form factor that range from 
to 3 (GeV/c)-i. 




FIG. 6: Helicity of the Dtt distribution for experimental events (points) and for MC simulation 
(histogram). The hatched distribution shows the background from the AE sideband region with 
a proper normalization, (a) corresponds to the D2 region \Mdt^ — 2.46| < O.lGeV/c^; (b) - Dq 
region {Mon - 2.30| < 0.1 GeV/c^. 

The helicities of the vrvr system in the M,r7r range of p, /2 and below the p where the 
broad resonance dominates are shown in Fig. [7| For the positive helicity range where the Dtt 
contribution is suppressed, a clear P-wave structure for p and D-wave for /2 is observed. The 
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branching ratios for the /o(600) is B(B^ foD^)B{fo tt+tt") = (1.75±0.26±0.35t|5i|) x 
10~^. This process can also have a contribution from some nonresonance background. The 
branching ratios for the p and the /2 cLFG clS follows: 

^ /2^°)i3(/2 ^ TT+TT-) = (1.10 ± 0.19 ± 0.21+°+^) X 10-^ 

Taking into account the PDG value of the branching fraction B{f2 tttt) = 0.848lo!oi3 ^^'^ 
Clebsch-Gordan coefficients,we obtain: 

B{B^ f2D°) = (1.95 ± 0.34 ± 0.38t°:^2) x 10-^ 

B{B° p^D^) = (2.91 ± 0.28 ± 0.331°;°^) x 10"^ 




FIG. 7: Helicity distribution for experimental events (points) and for MC simulation (histogram). 
The hatched distribution shows the background distribution from the AE sideband region with a 
proper normahzation. (a) corresponds to the p region |M^^ — 0.78| < 0.2GeV/c^; (b) - /2 region 
\M^^ - 1.20| < O.lGeV/c^; (c) - /o region M^^ < OMGeV/c^. 



D*^TT+TT- ANALYSIS. 

For D* reconstruction, the D*^ D^n^ decay is used with two decay modes of the 
£)0. jjo _^ K^Tc^ and D° — > K^n^Ti^'n-^ . The event distributions in AE and M^c are 
shown in Fig. |S| In each mode the number of signal events is obtained in a way similar to 
that described for the Dtttt selection. The observed signal yields of = 278 ± 23 and 
Nk3tt = 269 ± 29 for the Kn and Knirir modes, respectively, are consistent when efficiencies 
determined from MC, and the following D branching fractions are used: (3.80 ± 0.09)% for 
K~TT~^ and (7.46 ±0.31)% for T^'-tt+tt+tt-. The branching fraction of {D* —>■ Z}tt)tttt events, 
calculated from the weighted average of the values obtained for the two modes, is: 

B{B° D*\+n-) = (1.09 ± 0.08 ± 0.16) x 10'^, 

where the first error is statistical and the second is systematic. This measurement is about 
2.5(7 larger than the previous result Q of (0.62 ± 0.22 ± 0.22) x 10"^ We consider the 



previous measurement to be a statistical fluctuation. The systematic error contributions are 
listed in Table ITVl 

The background shape uncertainty and the effect of cut boundaries were estimated in the 
same way as for the Dtttt analysis. 
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FIG. 8: Distribution of D*^it^t: events in Ai? and Mbc plots with D reconstruction in K-k 

mode - (a), (c) and in Kdi-K mode - (b), (d). (a), (b) /S.E distributions, (c), (d) Mbc distributions. 



B — > D*TnT coherent amplitude analysis 



In this final state we have a decaying vector D* particle. There are two additional degrees 
of freedom and, in addition to the D*7i and tttt invariant masses squared (g^, qf), two other 
variables are needed to specify the final state. The variables are chosen to be the angle a 
between the pions from the D** and D* decay in the D* rest frame, and the azimuthal angle 
7 of the pion from the D* relative to the B D*!^ decay plane. For the case of the tttt 
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Source 


^ sys^ 


ijr[u , u ) 


O.O 


Tracking 


/I 


/I 


« 

u 


PID 


5 


MC 


3 


Background 


10 


Total 


14.7 



TABLE IV: Contributions to the systematic error for B D*inr. 

structure analysis another set of parameters can be chosen: vrvr square mass (qf)', the hehcity 
angle 6' of the vrvr-meson — the angle between the positive pion from this meson decay and 
the D* direction in the vrvr meson rest frame; the helicity angle a' of the D* meson — the 
angle between the pion from the D* decay and the vrvr-meson in the D* rest frame; and the 
angle 7' between the decay planes of the D* and the Trvr-meson. 

For further analysis, events satisfying the selection criteria described in the first section 
within the signal region ((AE + K(Mbc — MB))/aA£;)^ + ((Mbc — MB)/crMbJ^ < s are selected. 
The parameters a^E = llMeV/c^, ciMbc = 2.7 MeV/c^, k = 0.9 are obtained from a fit to 
experimental data, and the coefficient k takes into account the correlation between Mbc and 
AE. The parameter s is selected to be 4 for the mode with D — * Kn and 3 for D Kmiir 
to have a similar signal-to-background ratio. To understand the contribution and shape of 
the background, we use events in the sidebands {{AE ± 65 MeV + ^(Mbc - MB))/a^Ef + 
((Mbe - MB)/aMj^ < s. 

The D*7r and vrvr mass distributions for the signal and sideband events are shown in 
Fig. IHl and ^| There is a clear peak of the narrow states and Di with a negligible 
contribution of the broad states in the Dtt distribution. In the vrvr distribution the peaks 
of p and /2(1270) are clearly seen while the peak-like structure around 2.6 GeV/c^ is the 
reflection of the Dir angular distribution. In order to have the same Dalitz plot boundary 
for events from both the signal and the sideband regions as well as to decrease the smearing 
effect introduced by the detector resolution, mass-constrained fits of Dti to M£,*+ and Z}*vrvr 
to Mb are performed. 

To extract the amplitudes and phases for different intermediate states, an unbinned like- 
lihood fit in the four-dimensional phase space is performed. Assuming that the background 
distribution g^, a, 7) in the signal region has the same shape as in the AE sideband, 
we obtain the B{q^ , qf, a, 7) dependence from a fit of the sideband distribution to a smooth 
four- dimensional function. 

The number of background events in the signal region is normalized according to the 
relative areas of the signal and the sideband regions. The signal is parameterized as a sum 
of the amplitudes of an intermediate tensor {D2), two axial vector mesons {D[, Di), and 
three resonances p, /2 and /o(600) in the vrvr mode. For p and /2 there can be three different 
amplitudes depending on the relative polarizations of the decay products:74o, A± and A\\. 

Finally, the signal is expressed as follows: 

S{q\ ql a, 7) = ^2^(2) (g^ qf, a, 7) + a,e^'^^ A^'^\q^ ql a, 7) + a^e^'t'- A^^\q\ ql a, 7) 



16 




FIG. 9: Mass distribution of D*ir and tttt events, 
liistogram is the background distribution obtained 
simulation witli the amphtudes and parameters of 
fit. 



Points are the experimental data, the hatched 
from the sidebands, the open histogram is MC 
the intermediate resonances obtained from the 



+a,e^^^{{l - a|; - a'[)Ai'\q\ ql a, 7) + afje^^n ql a, 7) 

+a^^e^<A?(g^g?^a,7)) 

+aj,e^^f^{{l - af - g^^ a, 7) + al^e< Af\q\qla,^) 

+afle^^'^A'^t\q\qla,^)) 

+a},A^f"\q^,qla,^) + aps?. (4) 

The MC gives the resolution in invariant mass of about 1.9 MeV/c^, which is smaller than 
the resonance widths and can be neglected. 

Each resonance is described by a relativistic Breit-Wigner function with a width de- 
pending on q^. The angular dependence for each resonance corresponds to the spins of the 
intermediate and final state particles The amplitudes of 5 — > D* p and B D*f2 

decays have the following angular distributions: 

Aq ~ cos a' cos ^' 

y4|| ~ sin a' sin 6'' cos 7' 

~ sin a' sin ^' sin 7' (5) 

AI'' ~ cosa'(cos2 0' - 1/3) 
Aj^^ ~ sin a' sin 9' cos 9' cos 7' 

A-^f ~ sin a' sin 9' cos 9' sin 7' (6) 
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FIG. 10: The Dalitz plot of (a) signal events; (b) sideband events. 
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TABLE V: Comparison of the models with different resonances included in the fits. 

Table IVII shows the results of the fit for different models. If we remove the D[ meson 
from Eq. EJ the likelihood does not change significantly. If we remove the /2 instead, the 
likelihood increases by 41. Adding the phase space term does not improve the likelihood 
significantly. Replacing the /2 meson with an /o(1370) results in a worse likelihood value as 
shown in Table IVl 

The masses and widths of the vrvr resonances are fixed at their PDG values; M^.+ , and 
r^,+ are taken from the Dtttt fit; and Md[ = 2427MeV/c2, T%,^ = 384MeV/c2 have been 
taken from our measurement for _D**o The mass and width as well as the 

branching fractions and phases of the amplitudes were free parameters of the fit. 

In addition to the main minimum, there are several local minima that differ in —2 IuC/Cq 
by 2.1 to 25. The main minimum and local minima are treated in the same way as for Dtttt. 

Figs, ini and demonstrate the comparison of data and the MC simulation generated 
according to Eq. Q with parameters obtained from the fit. The distributions of the helicities 
of D** and D* and the angle 7 in different regions demonstrate reasonable agreement of 
the experimental data and fit result. 

For the Di meson we obtain the following parameters: 

= (2428.2 ± 2.9 ± 1.6 ± 0.6) MeV/c^ = (34.9 ± 6.61^;^ ± 4.1) MeV. 

These parameters are in good agreement with the CLEO measurement for D^: 
Mflo = (2425 ± 2 ± 2) MeV/c^, F^o = (26i? ± 4) MeV 0. 
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+41 


-1 


Br n*(10~'^) 


2 45± 42 


2 45± 42 


2 48± 43 


2 43± 41 


Brn, (10-^) 


0.908±0.145 
3.68± 0.60 


0.907±0.145 
3.71± 0.62 


0.837±0.139 
4.03± 0.84 


0.766±0.147 
3.63± 0.61 


Bra, (10-^) 


-0.197±0.584 
0.14± 0.13 




-0.316±0.670 
O.llib 0.12 


-0.121±0.556 
O.llib 0.14 


<l>ps 

5r„JlO-^) 


- 


- 


- 


2.594±0.551 
O.OOib 0.17 


BrJW'^) 

P \ ) 


2.566±0.333 
3.73± 0.87 


2.543±0.337 
3.78± 0.87 


2.032±0.326 
3.89± 0.96 


2.560±0.327 
3.74± 0.85 


Brf.ilO-^) 


0.440±0.413 
1.05± 0.37 


0.411±0.425 
1.05± 0.37 




0.429±0.453 
0.98± 0.35 


Brj.ilO-^) 


-2.263±0.646 
0.17± 0.11 


-2.190±0.643 
0.16± 0.11 


-2.823±0.498 
0.32± 0.17 


-2.181±0.597 
0.17± 0.11 


II 

< 


0.204±0.059 
0.067±0.038 
0.678±0.348 
0.730±0.058 
2.046±0.229 


0.198±0.058 
0.065±0.038 
0.686±0.351 
0.737±0.057 
2.031±0.229 


0.176±0.061 
0.105±0.042 
0.693±0.307 
0.719±0.059 
2.269±0.250 


0.211±0.059 
0.066±0.038 
0.624±0.358 
0.723±0.058 
2.030±0.225 


0^ 

«f 
0f 


0.623±0.137 
0.080±0.084 
-3.036±0.687 
0.297±0.137 
-0.895±0.489 


0.616±0.143 
0.092±0.091 
-2.983±0.672 
0.292±0.142 
-0.846±0.506 




0.646±0.142 
0.082±0.085 
-3.129±0.752 
0.273±0.143 
-0.926±0.551 



TABLE VI: The fit results for different sets of amplitudes. 



The preliminary results for the product of the branching ratios of the D**'s are the 
following: 

^ Din-) X B{Di D*\+) = (3.68 ± 0.60l™J^) x 10-^ 

^ D;+n-) X B{D;+ D*\+) = (2.45 ± 0.42;°i5;°:??) x 10-^ 
i3(fi" ^ D'+Ti-) X B{D'+ D*\+) = (0.14 ± 0.13 ± 0.12+°;?°) x 10"^ 
The last value is not statisticaly significant and corresponds to an upper limit: 

i3(5° ^ D'+Ti-) X B{D'+ D*\+) < 0.7 x 10"^ at 90% C.L. 

Including a contact term improves the likelihood but without high significance (see Table IVH) . 
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FIG. 11: The distribution of the data for the Di region: (a)-(c) {Md'tt - 2.41| < OMGeV/c^ and 
the region: \MD*n - 2.45| < 0.025 GeF/c^ (d)-(f). (a), (d) cosO - hehcity angle of D**; (b), 
(e) cosa - heUcity angle of D*; (c), (f) azimuthal angle 7. The points are experimental data, the 
histogram is MC with fitted parameters, and the hatched histogram is the background contribution 
(from sideband). 



The helicity of the vrvr system and the D* as well as the azimuthal angle 7' are plotted 
in Fig. 1121 for the M^,^ range of the p and the /2. 

The branching ratio for the /o(600) production is comparable with zero: B{B^ — *■ 
foD^)B{U vr+TT-) = (0.17 ± 0.11 ± 0.10t[j:J^) x 10^^ This contribution can also be 
regarded as some nonresonance background. 

The branching ratios observed for the p and the /2 are as follows: 

B{B° p^D*^) = (3.73 ± 0.87 ± 0.46l°;J^) x 10-^ 

B{B^ f2D*^)B{f2 TT+TT-) = (1.05 ± 0.37 ± 0.34+°;^^) x 10"^ 
Taking into account the branching fraction of B{f2 tt+tt^) we obtain: 

B{B^ f2D*^) = (1.86 ± 0.65 ± 0.60l°|°) x 10"^ 

In spite of the limited statistics, we can determine the contribution of different polariza- 
tion amplitudes for p and /2: 

ag = 0.73 ±0.06 ±0.10 ±0.09 
a(| = 0.20 ±0.06 ±0.03 ±0.10 
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FIG. 12: Helicity distribution for the p region: (a)-(c) \M^^ - .80| < 0.20 GeF/c^ and the 
region: \M^^ - 1.225| < 0.125 GeF/c^ (d)-(f). (a), (d) 008 6* - hehcity angle of D**; (b), (e) cos a 
- hehcity angle of D*; (c), (f) azimuthal angle 7. experimental events (points) and for fast MC 
simulation (histogram). The hatched distribution shows the background distribution from the AE 
sideband region with a proper normalization. 



Results and discussion 

The branching-fraction products obtained for the narrow {j = 3/2) resonances are similar 
to those obtained in the case of charged B decays as shown in Table IVlll The measured 
values of the branching fractions of the broad resonance and D'l' production in neutral 
B decay are, however, significantly lower than those for charged B decays. One possible 
explanation for this phenomenon is that for charged B decay to D**7t, the amplitude re- 
ceives contributions from both tree and the color suppressed diagrams as shown in Fig. ^1 
Production of D** via tree-diagrams is described by the Isgur-Wise functions T1/2 and T3/2. 
According to the sum rule |25l.l2&]. T1/2 ^ T3/2 and one would expect the suppression of the 
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Neutral B 


Charged B [17] 


B{B DiTT-)B{Di D*tt) 
B{B D^7r-)B{D^ Dtt) 


(3.68 ± 0.60t°:St°:^^) X 10-4 
(3.08 ± 0.33 ± 0.09l[!:J^) x 10"^ 


(-1 oinoinoinoN ^ i n-4 

(6.8 ± 0.7 ± 1.3 ±0.3) X 10-4 
(3.4 ±0.3 ±0.6 ±0.4) X 10-4 


B{B D[tt-)B{D[ D*tt) 
B{B D^tt)B{D*q Dtt) 


< 0.7 X 10-4 at 90% C.L. 

< 1.2 X 10-4 at 90% C.L. 


(5.0 ± 0.4 ± 1.0 ±0.4) X 10-4 
(6.1 ±0.6 ±0.9 ±1.6) X 10-4 



TABLE VIL Comparison of branching-fraction products for neutral and charged B decays. 




FIG. 13: Feynman diagrams for charged (a),(b) and neutral (c) B decays. 



broad state production. For the color suppressed diagrams, however, D**'s are produced by 
another mechanism and the amphtudes are characterized by the constants fD{3/2) and /d(i/2) 
and /d(3/2) /z)(i/2)- The production of the broad resonances Z)g° and D'^ in charged B 
decay is probably amplified by the color suppressed amplitude. 

CONCLUSION 

A study of neutral B to D^ti^ti^ and 0*^71^71^ decays has been presented. We have 
measured the total branching fractions of the three-body decays: B{B^ D^-k^-k^) = 
(1.07 ± 0.06 ± 0.10) X 10-3 and i3(S° D*^7i+7i-) = (1.09 ± 0.08 ± 0.16) x 10-^ and the 
two-body decay: i3(B0 ^ D*+7i-) = (2.30 ± 0.06 ± 0.19) x 10-^. 

The dynamics of these three-body decays has been studied. The D^ti^ti^ final state is 
described by the production of -D2^- with subsequent decays of D2 Dtx and Dp, D/2 
and a broad scalar (vrvr) structure. From a Dalitz distribution analysis we have obtained the 
branching fraction product for D2^: 

B{B^ D*+7i-) X B{D;+ D%+) = (3.08 ± 0.33 ± Omtoll) x 10"^ 

The values obtained for the mass and width of the tensor meson are: 

M^*+ = (2459.5 ± 2.3 ± 0.7tol)MeV /c^, F^.+ = (48.9 ± 5.4 ± 4.2 ± 1.9)MeV. 

The upper limit for the contribution of the scalar D^^ meson assuming its mass and 
width of Dq is: 

B{B° D*+Ti-) X B{D*+ D\+) < 1.2 x 10-4 at 90% C.L.. 
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The branching fractions for Dp and D/2 productions have been measured: 

^ p°L>°) = (2.91 ± 0.28 ± 0.331°:°^) x 10"^ 

B{B^ f2D^) = (1.95 ± 0.34 ± 0.381°;^^) ^ IQ-^. 

The last result represents a first observation. 

The D*7r7r final state is described by the production of DaTT and Ditt with subsequent 
decays of D** — > D*7r and D*p, D*f2. Prom a phase space analysis, we obtain the branching 
fractions product for D**^: 

5(5° ^ D+TT-) X B{D+ ^ L>*°7r+) = (3.68 ± 0.60l°:IJl°:^^) x 10-^ 
i3(5° ^ /^2*^7r-) X 5(^2*+ ^ D*°7r+) = (2.45 ± 0.42;0;^5+0;39) x 10"^ 
and set an upper limit on the production of the broad D[ resonance: 

B{B° D'+TT-) X B{D'+ D*%+) < 0.7 x 10"^ at 90% C.L.. 
For the Di meson mass and width we obtain the following values: 

= (2428.2 ± 2.9 ± 1.6 ± 0.6) MeV/c^ = (34.9 ± ± 4.1) MeV. 

The branching fraction of D*p and D*f2 has been measured: 

BiB'^ p'^D*'^) = (3.73 ± 0.87 ± 0.46l°:J^) x 10"^ 

B{B^ f2D*^) = (1.86 ± 0.65 ± 0.60t°|°) x 10-^ 

These are the first measurements of these processes. We also observe dominance of the 
longitudinal polarization amplitude for B — > D*p and B — > D*f2. 
All results are preliminary. 
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